Double-fc phase of the Dzyaloshinskii-Moriya helimagnet Ba 2 CuGe 2 07. 
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Neutron diffraction is used to re-investigate the magnetic phase diagram of the non- 
centrosymmetric tetragonal antiferromagnet Ba2CuGe207. A novel incommensurate double- k mag- 
netic phase is detected near the commensurate-incommensurate phase transition. This phase is 
stable only for field closely aligned with the 4- fold symmetry axis. The results emphasize the inad- 
equacy of existing theoretical models for this unique material, and points to additional terms in the 
Hamiltonian or lattice effects. 
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Dzyaloshinskii-Moriya (DM) antisymmetric exchange 
interactions [TJ [2] and their role in incommensurate mag- 
netic structures are once again at the forefront of con- 
densed matter research. One reason for the renewed 
interest is the discovery of novel topological skyrmion 
phases in the DM helimagnet MnSi and related com- 
pounds [SI [4] . The other is the role that DM interactions 
play in certain scenarios of the multiferroic effect [5j [6] . 
Arguably, no material has attracted more attention si- 
multaneously in both these contexts than the tetragonal 
DM helimagnet Ba2CuGe207. This compound is one of 
the very few systems where DM helimagnetism can co- 
exist with weak- ferromagnetic behavior [7]. Moreover, 
theory suggested that it can support stable skyrmions, 
similar to those in MnSi [7 . Ba2CuGe207 has an unusual 
soliton lattice structure in zero field |8] , and shows a very 
interesting behavior in applied fields [9-11 . While the 
material itself does not demonstrate spontaneous electric 
polarization, the isostructural and magnetically commen- 
surate Ba2CoGe207 [12] [13] is, in fact, an unconventional 
multiferroic material [13] . 

The most spectacular feature of Ba2CuGe207 is the 
commensurate-incommensurate (CI) transition that oc- 
curs in this material in applied magnetic fields [9] [TT] . 
Despite the uniqueness and fundamental importance of 
this compound, and years of systematic studies, its exact 
nature remains unresolved. Theory predicts a contin- 
uous transition with a divergent periodicity, that maps 
onto the classic Frenkel-Kontorova model [15]. In obvi- 
ous contradiction, experiments indicate a discontinuous 
reorientation of spins with a poorly defined "intermedi- 
ate state". The latter has been largely dismissed as a 
result of defect pinning [9] [TT], phase separation asso- 
ciated with a spin-flop- like discontinuous transition [16] 
or a poor alignment of the applied field relative to the 
unique 4- fold axis of the tetragonal structure [17 . In 
the present study we show that none of these explana- 
tions apply. In fact, the "intermediate state" is a novel 
double- k incommensurate phase that exists only for mag- 
netic fields almost perfectly aligned along the 4-fold axis. 
Being very unlike anything predicted for this material, 



the double-/c structure undermines all existing theoreti- 
cal descriptions of Ba2CuGe2 07. 

The crystal and magnetic structures of Ba2CuGe207 
are reviewed in detail in Ref.[TT] The magnetism is due to 
Cu 2+ ions that form a square lattice, cell-centered in the 
(a, b) plane of the tetragonal non-centric PA2\m struc- 
ture (Fig. lb in Ref. [11). Long-range helimagnetic or- 
dering occurs at Tn = 3.2 K. To a good approximation, 
the low-temperature spin arrangement is a cycloid with 
propagation vectors (1 ± £o> ±£o? 0), £o = 0.027, and the 
spins rotating in the (zbl, =pl, 0) planes, respectively, as il- 
lustrated in Fig. la of Ref. [TT]. (1, 0, 0) being the (ir, tt) 
point for the square lattice, nearest-neighbor spins are 
aligned almost antiparallel to each other. Due to tetrago- 
nal symmetry, one typically observes two equivalent mag- 
netic domains with propagation vectors (l + £o, £o> 0) and 
(1 + £o> —Co ? 0), respectively. Field-cooling in a magnetic 
field as small as a few tens of Gauss applied in the (a, b) 
plane ensures a single-domain structure [TO] . A more 
sizable in-plane field tends to re-orient the spin rotation 
plane of the entire cycloid to be perpendicular to the 
field direction. This is accompanied by a continuous and 
rather spectacular rotation of the propagation vector [10] . 
The above-mentioned CI transition occurs in magnetic 
fields applied along the c axis. It is preceded by a strong 
deformation of the planar cycloid, that in applied fields 
is better described as a "solition lattice" [9j [TT] . Beyond 
H c ~ 2.1 T the system is a commensurate AF with prop- 
agation vector (1, 0, 0). The focus of the present study is 
what happens just before this commensuration. 

A weak point of all previous neutron diffraction studies 
was a less than perfect alignment of the magnetic field 
with respect to the c axis. On the one hand, this mis- 
alignment potentially alters the nature of the phase tran- 
sition. On the other hand, measurements and the inter- 
pretation of the data are complicated by a continuous re- 
distribution of intensity between the two tetragonal do- 
mains. The latter is due to the non-zero in-plane compo- 
nent of the applied field, and may be history-dependent 
due to pinning. Thus, the first step in our strategy to 
clarify the issue of the "intermediate state" was to per- 



2 







1 10 100 800 








j^f I 








(b) 

(1,0£)_ 







0.5 1 1.5 2 2.5 3 

H(T) 

FIG. 1: (a) Neutron diffraction intensity measured in linear 
scans along the (l+£ cos(0— 7r/4), £ sin(0— 7r/4), 0) reciprocal- 
space direction for cj) = as a function of magnetic field ap- 
plied along the c axis, (b) Field dependence of the integrated 
intensity of the incommensurate satellite reflections at £ and 
3£ as well as of the commensurate magnetic Bragg peak at 
(1,0,0), obtained in Gaussian fits. For the measurements of 3£ 
the instrumental resolution was relaxed to optimize intensity. 

form a new set of diffraction measurements in a sample 
with an almost perfect alignment, while diligently follow- 
ing a field-cooling protocol for each field value. 

The data were taken on the TASP 3-axis spectrome- 
ter at Paul Scherrer Institut [18] . using neutrons of in- 
cident wavevector 1.3 A . The single crystal sample 
of mass 1.65 g was grown using the floating-zone tech- 
nique and mounted in a split-coil cryomagnet. A mis- 
alignment of less that 0.5° between the nominally vertical 
field and the crystallographic c axis at low temperatures 
was achieved. The data were taken after field-cooling 
from 6 K to 1.6 K, in sets of elastic scans along the 
(l + £cos(0 — 7r/4), £sin(0 — 7r/4), 0) directions in recipro- 
cal space, for <j> = and <p = tt/2. Here we have adopted 
the notation used in Ref. [10 , where <p is the angle be- 
tween the (1, 1,0) direction and the propagation vector. 
The background, primarily due to double scattering [9], 
was subtracted from all low-temperature scans. 

The result of these measurements are summarized in 
Figs.[l]and [2j With increasing magnetic field, we observe 
a gradual change of the magnetic period ( and a shift of 
intensity from the = domain to that with <j) = 7r/2. 
Even at H = 1.7 T a single-domain state is not achieved, 
though most of the intensity has shifted to <j) = tt/2. This 
confirms an almost-perfect alignment of our sample rela- 
tive to the c axis. A dramatic redistribution of intensity 
is seen at Hi = 1.9 T. At this point the incommensurate 
peaks in the <j) = tt/2 scans disappear. At the same time, 
those in <j) = scans abruptly gain intensity. In other 
words, at Hi the preferred incommensurate propagation 
vector abruptly rotates by tt/2. Moreover, an additional 
peak appears in all scans in the commensurate (1,0,0) 
position. The incommensurate peaks persist until it dis- 
appears at the CI transition at H c ~ 2.4 T. We identify 
the field range Hi < H < H c with the previously re- 
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FIG. 2: Neutron diffraction data as a function of magnetic 
field applied along the c axis for cj> — tt/2, analogous to Fig.[l] 
The intensity scale as given in Fig. [I] is used. 

ported intermediate state. 

The abrupt rotation of the incommensurate vector at 
Hi suggests that this state is not a phase-separated mix- 
ture of the commensurate and the low-field cycloidal 
phases, as previously speculated, but a distinct double-A: 
incommensurate structure. The same conclusion can be 
drawn from the behavior of higher-order harmonics, such 
as the ones at (1 + 3£ cos((/> — 7r/4), 3£ sin((/> — 7r/4), 0). At 
H < Hi, the 3rd harmonic increases with field (Fig. [| 
b), due to a distortion of the planar cycloid and the 
development of a soliton lattice [9 . In contrast, for 
Hi < H < H c the harmonics are absent. The incom- 
mensurate modulation is purely sinusoidal, without dis- 
tortion, and thereby distinct from that below Hi. 

A distinct incommensurate "cone" phase preceding to 
the CI transition was indeed predicted theoretically, for 
fields applied at an angle a > with respect to the 4- 
fold axis [17]. This phase is supposedly absent for a = 0, 
and its field range progressively increases with increas- 
ing a. To verify that the observed double- k structure is 
not due to a residual a < 0.5° field misalignment, we 
systematically investigated the behavior in canted fields. 

In the first such measurement the field was applied at 
an angle a — 5° relative to the c axis, in the (1,0,0) 
plane. Despite the much larger misalignment angle, we 
observe an almost unchanged behavior. These data are 
summarized in Fig. [3] a-c. Due to a sizable in-plane 
field component, only one domain, with </> ~ 0, survives 
beyond 0.5 T applied field. The same in-plane com- 
ponent induces a gradual rotation of the incommensu- 
rate peak around the (1,0,0) position, <p changing from 
(j) = at H = to - 15° at H = 1.9 T in agree- 
ment with Ref. [TO] . The direction of the propagation 
vector was determined by centering the incommensu- 
rate peaks at each field. The £-scans that make up the 
false-color plot in Fig. [3] b are thus performed along the 
(l+£cos((/> — 7r/4), £sin((/> — 7r/4), 0) direction with a field- 
dependent <j). As for zero tilt, at H = Hi ~ 1.9 T, the 
propagation vector abruptly rotates by precisely 90 de- 
grees {<j) — » 0+7r/2 ), and a commensurate component ap- 
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FIG. 3: Neutron diffraction data as a function of magnetic 
field for a = 5°. (a) Magnetic field dependence of 0. (b) Neu- 
tron diffraction intensity measured in linear scans along the 
(1 + £ cos(0 — 7r/4), £ sin((/> — 7r/4), 0) reciprocal-space direc- 
tion. The intensity scale as given in Fig. nl is used, (c) Field 
dependence of the integrated intensity of the incommensu- 
rate satellite reflections at £ as well as of the commensurate 
magnetic Bragg peak at (1,0,0), as obtained in gaussian fits. 

pears. The incommensurate peak vanishes at H c ~ 2.4 T. 

A key finding of this study is the double- k phase is 
not enhanced by a field misalignment with respect to the 
c axis, but is actually absent at large a angles. Data 
collected for a = 15° are shown in Fig.|4j Due to a much 
larger in-plane field component, the propagation vector 
rotates more rapidly and completely aligns itself with 
the crystallographic (1,0,0) direction: <p — >• 7r/4 for H > 
1 T. This gradual re-orientation is fully consistent with 
previous studies [10] . No abrupt re-orientation of the 
propagation vector is observed at any field. The behavior 
is continuous all the way to the CI transition at H c ~ 
2.5 T, with no sign of an additional phase at lower fields. 

In contrast to small-a geometries, where the CI tran- 
sition occurs in a sinusoidally modulated phase with no 
higher-order harmonics, for a = 15° we observe a strong 
distortion of sinusoidal modulation. Indeed, as clearly 
seen in Fig. [4j the third (3£), as well as the previously 
unobserved 2nd (2£) harmonics, drastically increase in 
intensity as the CI transition is approached. 

Note that in Fig. [4}}, at H > H c , the commensurate 
peak appears to be absent. This, however, is an instru- 
mental effect. It stands to reason that in the commen- 
surate AF phase, spins are aligned perpendicular to the 
field, i.e., along the (1,0,0) direction. Neutrons are only 
scattered by spin components that are transverse to the 
scattering vector. As a result, the (1,0,0) AF Bragg 
peak is extinguished by this polarization effect. Due to 
the restrictive geometry of a split-coil cryomagnet, mea- 
surements at equivalent wave vectors with more favorable 
polarization were not feasible. 

Preliminary data (not show) indicate virtually iden- 



FIG. 4: Neutron diffraction data as a function of magnetic 
field for a — 15°, analogous to Fig. [3] The intensity scale as 
given in Fig. [I] is used. 

tical behaviors for fields confined to the (1,1,0) plane 
[T9] . For a = 5° the propagation vector jumps by 90° 
at Hi ~ 1.9 T, while for a = 15° no reorientation is ob- 
served and the peak remain at <j) = all the way to the 
CI transition. We can therefore summarize all diffraction 
results as follows. In field applied almost parallel to the 4- 
fold axis, a novel sinusoidally modulated double- k phase 
appears just before the CI transition. When the field 
deviates from the high-symmetry direction, the double-Zc 
phase is absent. In this case the CI transition occurs in 
a complexly distorted single- k phase that continuously 
evolves from the zero field structure. 

Even the limited diffraction data allow us to make a 
good guess regarding the nature of the novel double- k 
phase. Since it has a commensurate antiferromagnetic 
component, the double- k structure is clearly not the con- 
ical phase discussed in Ref. [T?] . Instead, we propose 
that it can be approximated as an antiferromagnetic cone 
structure depicted in Fig. [5] The spins are mostly aligned 
in a commensurate AF pattern in the (a, b) plane. In 
addition, there is a small incommensurate precession of 
transverse spin components. We will assume that, due to 
tetragonal symmetry, this structure can freely rotate in 
the (a, b) plane, just like the soliton lattice [ID] . Similarly, 
we can expect that the orientation of the spin rotation 
plane is uniquely coupled to the propagation vector. 

In a field, any AF structure will always favor a spin-flop 
configuration. In particular, the dominant commensu- 
rate AF spin component of the proposed AF-cone struc- 
ture will tend to align itself perpendicular to any in-plane 
component of the applied field. As a result, the spin pre- 
cession plane will align parallel to the in-plane field com- 
ponent. The expected behavior for the planar soliton- 
lattice phase (Fig.^) is quite different. A small in-plane 
field will tend to align the spin rotation plane perpendic- 
ular to itself [ID] . Thus the observed abrupt 90° rotation 
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FIG. 5: Proposed magnetic structure of Ba2CuGe207for the 
almost AF cycloid, the soliton lattice and the AF-cone phase. 

of the propagation vector at Hi may signify a transition 
from the planar soliton-lattice phase at low fields to a 
AF-conical state at high fields. 

The proposed model for the double- k structure also 
gives a handwaving explanation of the observed absence 
of higher-order harmonics. Indeed, for the planar state 
realized below Hi , a distortion and formation of a soliton 
lattice is the only way it can respond to a field applied in 
the plane of spin rotation. In contrast, the AF cone state, 
can take advantage of the Zeeman energy by a canting of 
its dominant AF-commensurate component, leaving the 
small incommensurate modulation unperturbed. 

A further understanding of the observed double- k 
phase will require new theoretical input. Existing mod- 
els fail to predict it, and are clearly missing an impor- 
tant component of the physics of Ba2CuGe2 07. The 
new phase appears close to the CI transition, where the 
soliton lattice of the original model becomes infinitely 
soft. Therefore, the missing terms in the Hamiltonian 
need not be particularly strong to qualitatively alter the 
energy landscape. We can suggest two potential can- 
didates. First, many of previous theories focused only 
on the component of the DM vector D that lies in the 
tetragonal (a, b) plane. Due to crystallographic symme- 
try, this component is uniform from one Cu-Cu bond to 
the next, and is responsible for all the spiral structures 
in Ba2CuGe207. The component of D that is paral- 
lel to the four- fold axis (Fig. lb in Ref. [11]) has often 
been ignored. This component is sign-alternating and fa- 
vors a commensurate weak-ferromagnetic canting of all 
spins in one direction [20]. This term is known to be 
dominant in the similar K2V3O8 system [2T] , If active 
in Ba2CuGe207, it will distort any planar soliton-lattice 
phase, potentially giving rise to even-order harmonics, 
as those observed in this work for large a. Note that 
distortions due to magnetic field itself produce only odd 
harmonics [9]. We can speculate that this term could 
also play a role in stabilizing the double-A: structure. 

The other potential cause of the new phase is lattice 
effect. Indeed, all previous models were based on a con- 
tinuous mapping of the original Heisenberg spin Hamil- 
tonian. Near the CI transition though, such a description 
is bound to break down. In this theory, even as the dis- 
tance between solitons increases near H C1 the width of 
each soliton decreases. Near the critical field it will be- 



come comparable to the lattice spacing, at which point 
the theory will lose self-consistency. 

In summary, our results show that there is more to the 
physics of Ba2CuGe2 07 than previously thought. This 
circumstance may have a direct bearing on the stability 
of the predicted (but never observed) skyrmion lattice in 
this material [7j [20] , as well as the multiferroic effect. 
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